FRACTILE ANTENNA ARRAYS AND METHODS FOR PRODUCING A FRACTILE 

ANTENNA ARRAY 

FIELD OF THE INVENTION 

The present invention is directed to fractile antenna arrays and a method of 
5 producing a fractile antenna array with improved broadband performance. The present 
invention is also directed to methods for rapidly forming a radiation pattern of a fractile 
array. 

BACKGROUND OF THE INVENTION 

Fractal concepts were first introduced for use in antenna array theory by Kim and 

10 Jaggard. See, Y. Kim et al., "The Fractal Random Array, " Proc. IEEE, Vol. 74, No. 9, 
pp. 1278-1280, 1986. A design methodology was developed for quasi-random arrays 
based on properties of random fractals. In other words, random fractals were used to 
generate array configurations that are somewhere between completely ordered (i.e., 
periodic) and completely disordered (i.e., random). The main advantage of this technique 

15 is that it yields sparse arrays that possess relatively low sidelobes (a feature typically 
associated with periodic arrays but not random arrays) which are also robust (a feature 
typically associated with random arrays but not periodic arrays). More recently, the fact 
that deterministic fractal arrays can be generated recursively (i.e., via successive stages of 
growth starting from a simple generating array) has been exploited to develop rapid 

20 algorithms for use in efficient radiation pattern computations and adaptive beamforming, 
especially for arrays with multiple stages of growth that contain a relatively large number 
of elements. See, D. H. Werner et. al., "Fractal Antenna Engineering: The Theory and 
Design of Fractal Antenna Arrays, " IEEE Antennas and Propagation Magazine, Vol. 41, 



No. 5, pp. 37-59, October 1999. It was also demonstrated that fractal arrays generated in 
this recursive fashion are examples of deterministically thinned arrays. A more 
comprehensive overview of these and other topics related to the theory and design of 
fractal arrays may be found in D. H. Werner and R. Mittra, Frontiers in Electromagnetics 
5 (IEEE Press, 2000). 

Techniques based on simulated annealing and genetic algorithms have been 
investigated for optimization of thinned arrays. See, D. J. O'Neill, "Element Placement 
in Thinned Arrays Using Genetic Algorithms, " OCEANS '94, Oceans Engineering for 
Today's Technology and Tomorrows Preservation, Conference Proceedings, Vol. 2, pp. 

10 301-306, 199; G. P. Junker et al., "Genetic Algorithm Optimization of Antenna Arrays 
with Variable Interelement Spacings, " 1998 IEEE Antennas and Propagation Society 
International Symposium, AP-S Digest, Vol. 1, pp. 50-53, 1998; C. A. Meijer, 
"Simulated Annealing in the Design of Thinned Arrays Having Low Sidelobe Levels, " 
COMSIG'98, Proceedings of the 1998 South African Symposium on Communications 

15 and Signal Processing, pp. 361-366, 1998; A. Trucco et al., "Stochastic Optimization of 
Linear Sparse Arrays, " IEEE Journal of Oceanic Engineering, Vol. 24, No. 3, pp. 291- 
299, July 1999; R. L. Haupt, "Thinned Arrays Using Genetic Algorithms, " IEEE Trans. 
Antennas Propagat., Vol. 42, No. 7, pp. 993-999, July 1994. A typical scenario involves 
optimizing an array configuration to yield the lowest possible side lobe levels by starting 

20 with a fully populated uniformly spaced array and either removing certain elements or 
perturbing the existing element locations. Genetic algorithm techniques have been 
developed for evolving thinned aperiodic phased arrays with reduced grating lobes when 
steered over large scan angles. See, M. G. Bray et al., "Thinned Aperiodic Linear Phased 
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Array Optimization for Reduced Grating Lobes During Scanning with Input Impedance 
Bounds, " Proceedings of the 2001 IEEE Antennas and Propagation Society International 
Symposium, Boston, MA, Vol. 3, pp. 688-691, July 2001; M. G. Bray et al., "Matching 
Network Design Using Genetic Algorithms for Impedance Constrained Thinned Arrays, " 
5 Proceedings of the 2002 IEEE Antennas and Propagation Society International 
Symposium, San Antonio, TX, Vol. 1, pp. 528-531, June 2001; M. G. Bray et al, 
"Optimization of Thinned Aperiodic Linear Phased Arrays Using Genetic Algorithms to 
Reduce Grating Lobes During Scanning, " IEEE Transactions on Antennas and 
Propagation, Vol. 50, No. 12, pp. 1732-1742, Dec. 2002. The optimization procedures 

10 have proven to be extremely versatile and robust design tools. However, one of the main 
drawbacks in these cases is that the design process is not based on simple deterministic 
design rules and leads to arrays with non-uniformly spaced elements. 
SUMMARY OF THE INVENTION 

The present invention is directed to an antenna array, comprised of a fractile array 

15 having a plurality of antenna elements uniformly distributed along Peano-Gosper curve. 

The present invention is also directed to an antenna array comprised of an array 
having an irregular boundary contour. The irregular boundary contour comprises a plane 
tiled by a plurality of fractiles and the plurality of fractiles covers the plane without any 
gaps or overlaps. 

20 The present invention is also directed to a method for generating an antenna array 

having improved broadband performance. A plane is tiled with a plurality of non- 
uniform shaped unit cells of an antenna array. The non-uniform shape of the unit cells 
and the tiling of said unit cells are then optimized. 
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The present invention is also directed to a method for rapidly forming a radiation 
pattern of a fractile array. A pattern multiplication for fractile arrays is employed 
wherein a product formulation is derived for the radiation pattern of a fractile array for a 
desired stage of growth. The pattern multiplication for the fractile arrays is recursively 
5 applied to construct higher order fractile arrays. An antenna array is then formed based 
on the results of the recursive procedure. 

The present invention is also directed to a method for rapidly forming a radiation 
pattern of a Peano-Gosper fractile array. A pattern multiplication for fractile arrays is 
employed wherein a product formulation is derived for the radiation pattern of a fractile 
10 array for a desired stage of growth. The pattern multiplication for the fractile arrays is 
recursively applied to construct higher order fractile arrays. An antenna array is formed 
based on the results of the recursive procedure. 
BRIEF DESCRIPTION OF THE DRAWINGS 

The accompanying drawings, which are included to provide further understanding 
15 of the invention and are incorporated in and constitute part of this specification, illustrate 
embodiments of the invention and, together with the description, serve to explain the 
principles of the invention. 

In the drawings: 

Figures 1A-1C illustrate element locations and associated current distribution for 
20 stage 1, stage 2 and state 3 Peano-Gosper fractile arrays; 

Figures 2A-2C illustrate the first three stages in the construction of a self-avoiding 
Peano-Gosper curve; 
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Figures 3A-3C illustrate Gosper islands and their corresponding Peano-Gosper 
curves for (a) stage 1, (b) stage 2, and (c) stage 4; 

Figure 4 illustrates a plot of the normalized stage 3 Peano-Gosper fractile array 
factor versus for# for cp = 0°; 

Figure 5 illustrates a plot of the normalized stage 3 Peano-Gosper fractile array 
factor versus 9 for <p = 90°; 

Figure 6 illustrates a plot of the normalized stage 3 Peano-Gosper fractile array 
factor versus q> for 0 = 90° and d min = X ; 

Figure 7 illustrates a plot of the normalized stage 3 Peano-Gosper fractile array 
factor versus 9 for cp =26° and d min = X ; 

Figure 8 illustrates a plot of the normalized array factor versus 9 with (p= 0° for a 
uniformly excited 19x19 periodic square array; 

Figure 9 illustrates plots of the normalized array factor versus 9 with (p = 0° and 
d min = 2X for a stage 3 Peano-Gosper fractile array and a 19x19 square array; 

Figure 10 illustrates plots of the normalized array factor versus 9 for <p = 0° with 
main beam steered to 0 O = 45° and (p 0 = 0°; 

Figures 1 1 A-l 1C illustrate the structure of the Peano-Gosper fractile array based 
on tiling of Gosper islands; 

Figure 12 illustrates a graphical representation of a plane tiled with non-uniform 
shaped unit cells; 

Figure 13 is a flow chart illustrating a preferred embodiment of the invention; 
Figure 14 is a flow chart illustrating a preferred embodiment of the invention; and 
Figure 15 is a flow chart illustrating a preferred embodiment of the invention. 



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Figures 1 A-1C illustrate the antenna element locations and associated current 
amplitude excitations for a stage 1, stage 2 and stage 3 Peano-Gosper fracticle arrays\ 
where the antenna elements are distributed over a planar area (e.g., in free-space, over a 
5 geographical area, mounted on an Electromagnetic Band Gap (EBG) surface or an 
Artificial Magnetic Conducting (AMC) ground plane, mounted on an aircraft, mounted 
on a ship, mounted on a vehicle, etc.) A fractile array is defined as an array with a fractal 
boundary contour that tiles the plane without leaving any gaps or without overlapping, 
wherein the fractile array illustrates improved broadband characteristics. The numbers 1 

10 and 2 denote each antenna element's relative current amplitude excitation. The minimum 
spacing between antenna elements is assumed to be held fixed at a value of d min for each 
stage of growth. The antenna elements may be comprised of shapes and sizes of 
elements well know to those skilled in the art. Some examples of potential applications 
for this type of array are listed in Table 1. 

15 Table 1 



Application 


Frequency (GHz) 


Wavelength (cm) 


dmin (Cm) 


Broadband 
L - Band Array 


1-2 


30-15 


15 


Broadband 
S - Band Array 


2-4 


15-7.5 


7.5 


Broadband 

L-Band & S-Band Array 


1-4 


30-7.5 


7.5 


Broadband 
C - Band Array 


4-8 


7.5-3.75 


3.75 


Broadband 

S-Band & C-Band Array 


2-8 


15-3.75 


3.75 


Broadband 

X - Band Array 


8-12 


3.75-2.5 


2.5 


Broadband 

C-Band & X-Band Array 


4-16 


7.5-1.875 


1.875 


Broadband 


12-18 


2.5-1.667 


1.667 
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K u - Band Array 








Broadband 
K - Band Array 


18-27 


1.667- 1.111 


1.111 


Broadband 

K a - Band Array 


27-40 


1.111-0.75 


0.75 


Broadband 

K u -, K-, & K a -Band Array 


12-48 


2.5-0.625 


0.625 


Broadband Millimeter 
Wave Array 


40-160 


0.75-0.1875 


0.1875 



Referring to Figures 2A-2C, the first three stages in the construction of a Peano- 
Gosper curve are illustrated. The generator at stage P=l, Figure 2 A, is first scaled by the 
appropriate expansion factor 5 to obtain the stage P-2 (Figure2B) construction of the 

5 Peano-Gosper curve. The expansion factor 8 is defined in equation 13, below, for a 
Peano-Gosper array. The next step in the construction process is to then replace each of 
the seven segments of the scaled generator by an exact copy of the original generator 
translated and rotated as shown in Figure 2B. This iterative process may be repeated to 
generate Peano-Gosper curves up to an arbitrary stage of growth P. Figures 3A-3C show 

10 stage 1, stage 2, and stage 4 Gosper islands bounding the associated Peano-Gosper curves 
which fill the interior. 

Higher-order Peano-Gosper fractal arrays (i.e., arrays with P>1) are recursively 
constructed using a formula for copying, scaling, rotating, and translating of the 
generating array defined at stage 1 Equations 1-14, below, are used for this 

15 recursive construction procedure. Figures 1 A-1C illustrate a graphical representation of 
the procedure. The array factor (i.e., radiation pattern) for a stage P Peano-Gosper fractal 
array is expressed in terms of the product of P 3x3 matrices which are pre-multiplied by a 
vector A and post-multiplied by a vector C. 
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AF P (Q,<p) = AB P C 



(1) 



where 
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A = [a, a 2 fl 3 ] 



a, = 2 cos 



* vv * mm 



sin0 cos($? - + (/* - 



27T 



C = 



^> = [/// 1(3x3) 





arctan 


y« 

K X n ) 


r„=< 




= 0 






( \ 

y n 

\ X n) 




arctan 


Ik 

<Pj = U-» t 



>*„>o 



+ n , x n < 0 



(2) 
(3) 

(4) 
(5) 

(6) 

(7) 
(8) 

(9) 



(10) 



(11) 
(12) 
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8 = 



k = 



2 sin a 
In 



(13) 
(14) 



where X is the free-space wavelength of the electromagnetic radiation produced by the 
fractile array. The selection of constants and coefficients are within the ordinary skill of 
5 the art. The values of N {j required in (8) are found from 



n=K1 = 



{1,3,5,6} {2} {4,7} 
{4,7} {1,3,5,6} {2} 
{2} {4,7} {1,3,5,6} 



(15) 



Expressions for (x n ,y n ) in terms of the array parameters d min ,a , and 8 for n = 1 - 7 are 
listed in Table 2. 

Table 2 



n 


X n 




1 


0.5d m i n (cosa-5) 


-0.5d m j n sina 


2 


0 


0 


3 


d m in(0.58-1.5cosa) 


1.5d m j n sina 


4 


d m in(0.58-2cosa-0.5cos(n/3+a)) 


d m in(0.5sin(ji/3+a)+2sina) 


5 


dmin(0.58- 1 .5cosa-cos(7i/3+a)) 


d m in(sin(7i/3+a)+l .5sina) 


6 


d m in(0.58-0.5cosa-cos(7i/3+a)) 


d m i n (sin(7i/3+a)+0.5sina) 


7 


d min (0.58-0.5cos(7i/3+a)) 


0.5d m inSin(7i/3+a) 
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With reference to Figure 4, a plot of the normalized array factor versus 9 for a 
stage 3 Peano-Gosper fractal array with <p = 0°is illustrated. Curve 410 represents the 
corresponding radiation pattern slices for the Peano-Gosper array with element spacings 
of d min = X . Curve 420 represents radiation pattern slices for a Peano-Gosper array with 
5 element spacings of d min =X/2. Likewise with reference to Figure 5, a plot of the 

normalized array factor versus 8 for a stage 3 Peano-Gosper fractal array with cp = 90° is 
illustrated. Curve 510 represents the corresponding radiation pattern slices for the Peano- 
Gosper array with element spacings of </ min = X and curve 520 represents radiation 
pattern slices for a Peano-Gosper array with element spacings of d min =X/2 .For Figures 
10 4 and 5, the angle cp is measured from the x-axis and the angle 0 is measured from the z- 
axis. 

With reference to Figure 6, a plot of the normalized array factor versus (p for a 
stage 3 Peano-Gosper fractile array where d mln = X, 0 = 90°, and 0° <<p< 360°. Figure 
6 demonstrates the absence of grating lobes present anywhere in the azimuthal plane of 

15 the Peano-Gosper fractile array, even with antenna elements spaced one-wavelength 
apart. The plot shows that the highest sidelobes in the azimuthal plane are 23.85 dB 
down from the main beam at 9 = 0°. The plot shown in Figure 6 also indicates that one 
of these sidelobes is located at the point corresponding to 6 = 90° and <p - 26°. A plot of 
the normalized array factor versus 0 for this Peano-Gosper factile array with (p = 26° and 

20 d m[n = X is shown in Figure 7. 

The plots illustrated in Figures 6 and 7 demonstrate that, for Peano-Gosper 
fractile arrays, no grating lobes appear in the radiation pattern when the minimum 



element spacing is changed from a half-wavelength to at least a full-wavelength. This 
results from the arrangement (i.e., tiling) of parallelogram cells in the plane forming an 
irregular boundary contour by filling a closed Koch curve. 

This result is in contrast to a uniformly excited periodic 19x19 square array, of 
5 comparable size to the stage 3 Peano-Gosper fractile array, containing a total of 344 
antenna elements. Referring to Figure 8, plots of the normalized array factor versus 
0and (p= 0° for the 19x19 periodic square array are illustrated for antenna element 
spacings of d m - n = d = X 1 2 , curve 820, and d min = d = X , curve 810 where the main 
beam orientation is 0 o = 0°and <p Q = 0°. A grating lobe is clearly visible for the case in 

10 which the elements are periodically spaced one wavelength apart. 

Referring to Figure 9, a plot 910 of the stage 3 Peano-Gosper fractile array factor 
versus 0 with (p= 0° is illustrated for the case where the minimum spacing between 
antenna elements is increased to two wavelengths (i.e., d min = 2X ). In contrast, a plot 
920 of the array factor versus 9 with q>= 0° for a uniformly excited 19x19 square array 

1 5 with elements spaced two wavelengths apart is also illustrated. Two grating lobes are 
clearly identifiable in the radiation pattern of the conventional 19x19 square array. 

The maximum directivity of a Peano-Gosper fractile array differs from that of a 
convention 19x19 square array. This value is calculated by expressing the array factor 
for a stage P Peano-Gosper fractile array with N p elements in an alternative form given 

20 by: 

^(^) = Z 7 * ex PW (16) 

n=l /7=1 
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where I n and /?„ represents the excitation current amplitude and phase of the n th element 
respectively,/; is the horizontal position vector for the n th element with magnitude r n and 

angle (p Hy and n is the unit vector in the direction of the far-field observation point. An 
expression for the maximum directivity of a broadside stage P Peano-Gosper fractile 
array, where the main bean is directed normal to the surface of the planar array, of 
isotropic sources may be readily obtained by setting fi n = 0 in (16) and substituting the 
result into 



|2 



D P = ' - ^ Lax (17) 



1 4.71 n 

— l\\AF p (0,<p)\ 2 sm6d6d<p 



*> K 0 0 

This leads to the following expression for the maximum directivity given by: 

(N p V 

E'. 

n=\ m~2 n = \ 

where 

^ = — J jcos(fc|r m - Fjsinflcos^-^)) smMpdO (19) 

Z7r 0 0 

and 0> m „ represents the polar angle measured from the x-axis to the vector r m = r m - r n . 
The inner integral in (19) may be shown to have a solution of the form 
1 2 " 

— Jcos(*|? m - ?Jsin0cos(p-flj)cty> = J 0 (k\r m - rjsinfl) (20) 
Substituting (20) into (19) yields 
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0 

The following integral relation (22) is then introduced 

K_ 

\ 

)J o (xsm0)sm0d0 = - 



(21) 



sin£ 

x 



(22) 



which may be used to show that (21) reduces to 
5 5=2- 



k K~r n \ 

Finally, substituting (23) into (18) results in 



(23) 







2 


IA 


m=2n=l 









(24) 



Table 3 includes the values of maximum directivity, calculated using (24), for 
several Peano-Gosper fractile arrays with different minimum element spacings d^ and 

10 stages of growth P. Table 4, furthermore, provides a comparison between the maximum 
directivity of a stage 3 Peano-Gosper fractile array and that of a conventional uniformly 
excited 19x19 planar square array. These directivity comparisons are made for three 
different values of antenna element spacings (i.e., d min = A/4, </ min = XI2 9 and d mln = X). 
Where the element spacing is assumed to be d min = XI 4 and d min = XI 2 , the maximum 

15 directivity of the stage 3 Peano-Gosper fractile array and the 19x19 square array are 
comparable. However, when the antenna element spacing is increased to d min = X , the 
maximum directivity for the stage 3 Peano-Gosper fractile array is about 10 dB higher 

Table 3 
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Minimum Spacing 


Stage Number 


Maximum Directivity 




P 


D p (dB) 


0.25 


1 


3.58 


0.25 


2 


12.15 


0.25 


3 


20.67 


0.5 


1 


9.58 


0.5 


2 


17.90 


0.5 


3 


26.54 


1.0 


1 


9.52 


1.0 


2 


21.64 


1.0 


3 


31.25 



than the 19x19 square array. This is because the maximum directivity for the stage 3 
Peano-Gosper fractile array increases from 26.54 dB to 31.25 dB when the antenna 
element spacing is changed from a half-wavelength to one-wavelength respectively. In 
5 contrast, the maximum directivity for the 19x19 square array drops from 27.36 dB down 
to 21.27 dB. The drop in value of maximum directivity for the 19x19 square array may 
result from the appearance of grating lobes in the radiation pattern. 



Table 4 



Element Spacing 
dmin / ^ 


Maximum Directivity (dB) 


Stage 3 Peano-Gosper Array 


19x19 Square Array 


0.25 


20.67 


21.42 


0.5 


26.54 


27.36 


1.0 


31.25 


21.27 
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Referring to Figure 10, a plot of the normalized array factor versus 0 for q> = 0° is 
illustrated where the main beam of the Peano-Gosper fractal array is steered in the 
5 direction corresponding to 6 0 = 45 °and q> 0 = 0°. The antenna element phases for the 
Peano-Gosper fractal array are chosen according to 

Pn = -K s in Go cos (^ " <P n ) (25) 
Curve 1010 shows the normalized array factor for a stage 3 Peano-Gosper fractal array 
where the minimum spacing between elements is a half-wavelength and curve 1020 
10 shows the normalized array factor for a conventional 19x19 uniformly excited square 
array with half-wavelength element spacings. This comparison demonstrates that the 
Peano-Gosper fractile array is superior to the 19x19 square array in terms of its overall 
sidelobe characteristics in that more energy is radiated by the main bean rather than in 
undesirable directions. 

15 Referring to Figures 1 1 A-l 1C, Peano-Gosper arrays are self-similar since they 

may be formed in an iterative fashion such that the array at stage P is composed of seven 
identical stage P-l sub-arrays (i.e., they consist of arrays of arrays). For example in 
Figure 1 IB, the stage 3 Peano-Gosper array is composed of seven stage 1 sub-arrays, 
Figure 1 1 A. Likewise, the stage 4 Peano-Gosper array, Figure 1 1C, consists of seven 

20 stage 2 sub-arrays, and so on. This arrangement of sub-arrays through an iterative 
process lends itself to a convenient modular architecture whereby each of these sub- 
arrays may be designed to support simultaneous multibeam and multifrequency 
operation. 
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This invention also provides for an efficient iterative procedure for calculating the 
radiation patterns of these Peano-Gosper fractal arrays to arbitrary stage of growth P 
using the compact product representation given in equation (6). This property may be 
useful for applications involving array signal processing. This procedure may also be 
5 used in the development of rapid (signal processing) algorithms for smart antenna 
systems. 

With reference to Figure 12, a graphical representation of a plane tiled with non- 
uniform shaped unit cells is illustrated. This invention also provides for a method of 
generating any planar or conformal array configuration that has an irregular boundary 
10 contour and is composed of unit cells (i.e., tiles) having different shapes. With reference 
to Figure 13, a flow chart is shown illustrating a method of the present invention for 
generating an antenna array having improved broadband performance wherein the 
antenna array has an irregular boundary contour. In step 13 10, a plane is tiled with a 
plurality of non-uniform shaped unit cells of an antenna array. In step 1320, the non- 
15 uniform shape of the unit cells are optimized. In step 1330, the tiling of said unit cells are 
optimized. The optimization may be performed using genetic algorithms, particle swarm 
optimization or any other type of optimization technique. 

With reference to Figure 14, a flow chart is shown illustrating a method of the 
present invention for rapid radiation pattern formation of a fractile array. In step 1410, a 
20 factile array initiator and generator are provided. In step 1420, the generator is 

recursively applied to construct higher order fractile arrays. In step 1430, a fractile array 
is formed based on the results of the recursive procedure. 
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With reference to Figure 15, a flow chart is shown illustrating a method of the 
present invention for rapid radiation pattern formation of a Peano-Gosper fractile array. 
In step 1510, a pattern multiplication for fractile arrays is employed wherein a product 
formulation for the radiation pattern of a fractile array for a desired stage of growth is 
5 derived. In step 1520, the pattern multiplication procedure is recursively applied to 
construct higher order fractile arrays. In step 1530, an antenna array is formed based on 
the results of the recursive procedure. 

The present invention may be embodied in other specific forms without departing 
from the spirit or essential attributes of the invention. Accordingly, reference should be 
10 made to the appended claims, rather than the foregoing specification, as indicating the 
scope of the invention. Although the foregoing description is directed to the preferred 
embodiments of the invention, it is noted that other variations and modification will be 
apparent to those skilled in the art, and may be made without departing from the spirit or 
scope of the invention. 
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